Non-random mating within an Island rookery of Hawaiian hawksbill turtles: demographic discontinuity at a small coastline scale.

Hawksbill sea turtles (Eretmochelys imbricata) from the Hawaiian archipelago form a small and genetically isolated population, consisting of only a few tens of individuals breeding annually. Most females nest on the island of Hawai'i, but little is known about the demographics of this rookery. This study used genetic relatedness, inferred from 135 microhaplotype markers, to determine breeding sex-ratios, estimate female nesting frequency and assess relationships between individuals nesting on different beaches. Samples were collected during the 2017 nesting season and final data included 13 nesting females and 1002 unhatched embryos, salvaged from 41 nests, of which 13 had no observed mother. Results show that most females used a single nesting beach laying 1-5 nests each. From female and offspring alleles, the paternal genotypes of 12 breeding males were reconstructed and many showed high relatedness to their mates. Pairwise relatedness of offspring revealed one instance of polygyny but otherwise suggested a 1 : 1 breeding-sex ratio. Relatedness analysis and spatial-autocorrelation of genotypes indicate that turtles from different nesting areas do not regularly interbreed, suggesting that strong natal homing tendencies in both sexes result in non-random mating across the study area. Complexes of nearby nesting beaches also showed unique patterns of inbreeding across loci, further indicating that Hawaiian hawksbill turtles have demographically discontinuous nesting populations separated by only tens of km.

Resolving fine-scale population structure and fishery exploitation using sequenced microsatellites in a northern fish.

The resiliency of populations and species to environmental change is dependent on the maintenance of genetic diversity, and as such, quantifying diversity is central to combating ongoing widespread reductions in biodiversity. With the advent of next-generation sequencing, several methods now exist for resolving fine-scale population structure, but the comparative performance of these methods for genetic assignment has rarely been tested. Here, we evaluate the performance of sequenced microsatellites and a single nucleotide polymorphism (SNP) array to resolve fine-scale population structure in a critically important salmonid in north eastern Canada, Arctic Charr (Salvelinus alpinus). We also assess the utility of sequenced microsatellites for fisheries applications by quantifying the spatial scales of movement and exploitation through genetic assignment of fishery samples to rivers of origin and comparing these results with a 29-year tagging dataset. Self-assignment and simulation-based analyses of 111 genome-wide microsatellite loci and 500 informative SNPs from 28 populations of Arctic Charr in north-eastern Canada identified largely river-specific genetic structure. Despite large differences (~4X) in the number of loci surveyed between panels, mean self-assignment accuracy was similar with the microsatellite loci and the SNP panel (>90%). Subsequent analysis of 996 fishery-collected samples using the microsatellite panel revealed that larger rivers contribute greater numbers of individuals to the fishery and that coastal fisheries largely exploit individuals originating from nearby rivers, corroborating results from traditional tagging experiments. Our results demonstrate the efficacy of sequence-based microsatellite genotyping to advance understanding of fine-scale population structure and harvest composition in northern and understudied species.

Chromosome polymorphisms track trans-Atlantic divergence and secondary contact in Atlantic salmon.

Pleistocene glaciations drove repeated range contractions and expansions shaping contemporary intraspecific diversity. Atlantic salmon (Salmo salar) in the western and eastern Atlantic diverged >600,000 years before present, with the two lineages isolated in different southern refugia during glacial maxima, driving trans-Atlantic genomic and karyotypic divergence. Here, we investigate the genomic consequences of glacial isolation and trans-Atlantic secondary contact using 108,870 single nucleotide polymorphisms genotyped in 80 North American and European populations. Throughout North America, we identified extensive interindividual variation and discrete linkage blocks within and between chromosomes with known trans-Atlantic differences in rearrangements: Ssa01/Ssa23 translocation and Ssa08/Ssa29 fusion. Spatial genetic analyses suggest independence of rearrangements, with Ssa01/Ssa23 showing high European introgression (>50%) in northern populations indicative of post-glacial trans-Atlantic secondary contact, contrasting with low European ancestry genome-wide (3%). Ssa08/Ssa29 showed greater intrapopulation diversity, suggesting a derived chromosome fusion polymorphism that evolved within North America. Evidence of potential selection on both genomic regions suggests that the adaptive role of rearrangements warrants further investigation in Atlantic salmon. Our study highlights how Pleistocene glaciations can influence large-scale intraspecific variation in genomic architecture of northern species.

Environmental extremes drive population structure at the northern range limit of Atlantic salmon in North America.

Conservation of exploited species requires an understanding of both genetic diversity and the dominant structuring forces, particularly near range limits, where climatic variation can drive rapid expansions or contractions of geographic range. Here, we examine population structure and landscape associations in Atlantic salmon (Salmo salar) across a heterogeneous landscape near the northern range limit in Labrador, Canada. Analysis of two amplicon-based data sets containing 101 microsatellites and 376 single nucleotide polymorphisms (SNPs) from 35 locations revealed clear differentiation between populations spawning in rivers flowing into a large marine embayment (Lake Melville) compared to coastal populations. The mechanisms influencing the differentiation of embayment populations were investigated using both multivariate and machine-learning landscape genetic approaches. We identified temperature as the strongest correlate with genetic structure, particularly warm temperature extremes and wider annual temperature ranges. The genomic basis of this divergence was further explored using a subset of locations (n = 17) and a 220K SNP array. SNPs associated with spatial structuring and temperature mapped to a diverse set of genes and molecular pathways, including regulation of gene expression, immune response, and cell development and differentiation. The results spanning molecular marker types and both novel and established methods clearly show climate-associated, fine-scale population structure across an environmental gradient in Atlantic salmon near its range limit in North America, highlighting valuable approaches for predicting population responses to climate change and managing species sustainability.

Annotated mitochondrial genome assemblies for two sand lances (genus: Ammodytes) from the northwest Atlantic.

Complete mitochondrial genomes of two northwest Atlantic sand lances (Ammodytes americanus and Ammodytes dubius) were sequenced, assembled, and annotated. Both genomes were 16 519 bp in length and were differentiated by a genetic distance of only 0.01. Furthermore, mitochondrial gene annotations were identical for both species. Phylogenetic analysis revealed that divergence between the two species was shallow, relative to other members of the genus.

Complete mitochondrial genomes for Icelus spatula, Aspidophoroides olrikii and Leptoclinus maculatus: pan-Arctic marine fishes from Canadian waters.

Three Arctic marine fishes Icelus spatula, Aspidophoroides olrikii and Leptoclinus maculatus have been identified as target species for investigating the effects of ocean warming on population patterns in high-latitude marine habitats around Canada. In preparation for this research, we have resolved whole mitochondrial genome sequences of 16 384, 17 200 and 16 384 bp for each species, respectively. GC content for each species was 47.5%, 44.2% and 45.3%, respectively. Mitogenome gene composition included 13 protein-encoding genes, 2 rRNA and 22 tRNA genes, for I. spatula and L. maculatus, consistent with other teleosts. Only 20 tRNA genes were annotated for A. olrikii, because tRNA-Pro and tRNA-Thr are poorly characterized and aberrantly located in this species.

Limited contemporary gene flow and high self-replenishment drives peripheral isolation in an endemic coral reef fish.

Extensive ongoing degradation of coral reef habitats worldwide has lead to declines in abundance of coral reef fishes and local extinction of some species. Those most vulnerable are ecological specialists and endemic species. Determining connectivity between locations is vital to understanding recovery and long-term persistence of these species following local extinction. This study explored population connectivity in the ecologically-specialized endemic three-striped butterflyfish (Chaetodon tricinctus) using mt and msatDNA (nuclear microsatellites) to distinguish evolutionary versus contemporary gene flow, estimate self-replenishment and measure genetic diversity among locations at the remote Australian offshore coral reefs of Middleton Reef (MR), Elizabeth Reef (ER), Lord Howe Island (LHI), and Norfolk Island (NI). Mt and msatDNA suggested genetic differentiation of the most peripheral location (NI) from the remaining three locations (MR, ER, LHI). Despite high levels of mtDNA gene flow, there is limited msatDNA gene flow with evidence of high levels of self-replenishment (≥76%) at all four locations. Taken together, this suggests prolonged population recovery times following population declines. The peripheral population (NI) is most vulnerable to local extinction due to its relative isolation, extreme levels of self-replenishment (95%), and low contemporary abundance.

Observations of migrant exchange and mixing in a coral reef fish metapopulation link scales of marine population connectivity.

Much progress has been made toward understanding marine metapopulation dynamics, largely because of multilocus microsatellite surveys able to connect related individuals within the metapopulation. However, most studies are focused on small spatial scales, tens of kilometers, while demographic exchange at larger spatial scales remains poorly documented. Additionally, many small-scale demographic studies conflict with broad-scale phylogeographic patterns concerning levels of marine population connectivity, highlighting a need for data on more intermediate scales. Here, we investigated demographic recruitment processes of a commercially important coral reef fish, the bluespine unicornfish (Naso unicornis) using a suite of mitochondrial DNA (mtDNA) and microsatellite markers. Sampling for this study ranged across the southern Marianas Islands, a linear distance of 250 km and included 386 newly settled postlarval recruits. In contrast with other studies, we report that cohorts of recruits were genetically homogeneous in space and time, with no evidence of temporally stochastic sweepstakes reproduction. The genetic diversity of recruits was high and commensurate with that of the adult population. In addition, there is substantial evidence that 2 recruits, separated by 250 km, were full siblings. This is the largest direct observation of dispersal to date for a coral reef fish. All indications suggest that subpopulations of N. unicornis experience high levels of demographic migrant exchange and metapopulation mixing on a spatial scale of hundreds of kilometers, consistent with high levels of broad-scale genetic connectivity previously reported in this species.

Limited ecological population connectivity suggests low demands on self-recruitment in a tropical inshore marine fish (Eleutheronema tetradactylum: Polynemidae).

The diversity of geographic scales at which marine organisms display genetic variation mirrors the biophysical and ecological complexity of dispersal by pelagic larvae. Yet little is known about the effect of larval ecology on genetic population patterns, partly because detailed data of larval ecology do not yet exist for most taxa. One species for which this data is available is Eleutheronema tetradactylum, a tropical Indo-West Pacific shorefish. Here, we use a partial sequence mitochondrial cytochrome oxidase subunit 1 (COI) marker and five microsatellite loci to survey the genetic structure of E. tetradactylum across northern Australia. Structure was found throughout the range and isolation by distance was strong, explaining approximately 87 and 64% of the genetic variation in microsatellites and mtDNA, respectively. Populations separated by as little as 15 km also showed significant genetic structure, implying that local populations are mainly insular and self-seeding on an ecological time frame. Because the larvae of E. tetradactylum have lower swimming performance and poor orientation compared with other tropical fishes, even modest larval abilities may permit self-recruitment rather than passive dispersal.

High population connectivity across the Indo-Pacific: Congruent lack of phylogeographic structure in three reef fish congeners.

We used the mitochondrial control region and a comparative approach to study the genetic population structure of two surgeonfishes, Naso brevirostris and Naso unicornis, across their Indo-central Pacific ranges. Our purpose was to compare our results with those of a previous study of Naso vlamingii [Klanten, S.O., van Herwerden, L., Choat J.H., 2007. Extreme genetic diversity and temporal rather than spatial partitioning in a widely distributed coral reef fish. Mar. Biol. 150, 659-670] another widely distributed Indo-central Pacific Naso species. We found no evidence of a barrier to gene flow between the Indian and Pacific Oceans for either species, consistent with what was shown for N. vlamingii. Overall, both target species lacked spatial population partitions and probably have complex patterns of gene flow on several spatial scales. Despite the lack of geographic population structure distinct clades were observed in N. brevirostris, similar to those found in N. vlamingii. Coalescence times for intraspecific clades of N. brevirostris and N. vlamingii approximate each other, suggesting parallel evolutionary histories. A bimodal mismatch distribution in N. brevirostris indicates that a biogeographic barrier separated N. brevirostris populations sometime during its species history. Naso unicornis, in contrast, lacked genetic structure of any kind, although it has what could represent a single surviving clade. Congruent lack of spatial population structure among all three species suggest that such patterns are not due to stochastic processes of DNA mutation and are most likely driven by ecological and environmental factors.